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Abstract 

Analytical solutions for the time-dependent autocorrelation function of the classical and 
quantum mechanical spin dimer with arbitrary spin are presented and compared. For large 
spin quantum numbers or high temperature the classical and the quantum dimer become 
more and more similar, yet with the major difference that the quantum autocorrelation 
function is periodic in time whereas the classical is not. 
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1 Introduction and summary 



There is a growing interest in the magnetic properties of synthesized molecules 
[1-4] containing relatively small numbers of paramagnetic ions. With the ability 
to control the placement of magnetic moments of diverse species within stable 
molecular structures, one can test basic theories of magnetism and even begin to 
explore the design of novel systems that offer the prospect of useful applications. 
Most species of organic -based molecular magnets exhibit very weak intermolecular 
magnetic interactions, so that measurements performed on a bulk sample actually 
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reflect intramolecular interactions only. The magnetic interaction appears to be well 
described by the Heisenberg model with isotropic, nearest-neighbor exchange. A 
key quantity is the time- and temperature-dependent correlation function for pairs 
of magnetic moments, as it serves as the basic ingredient for understanding diverse 
dynamical phenomena, such as inelastic neutron scattering [5] and spin lattice re- 
laxation [6]. 

The present study is motivated by a desire to achieve a deeper understanding of spin 
dynamics in the Heisenberg model, especially concerning the trends that occur in 
arrays of N interacting moments (individual spins s) for increasing values of both 
N and s. The classical Heisenberg model turns out to provide accurate quantita- 
tive results for static properties, such as magnetic susceptibility, down to thermal 
energies of the order of the exchange coupling [7,8]. It is quite easy to establish 
the connection of that model, for static properties, to the corresponding quantum 
model for arbitrary s. However, considerable care is required to successfully link 
up with classical Heisenberg spin dynamics starting from quantum Heisenberg spin 
dynamics. 

In this article we present the analytical form of the time-dependent equilibrium 
autocorrelation function of the quantum mechanical dimer with general spin s. The 
trends for increasing s are explored in some detail and in particular we compare 
the quantum results with the exact analytical result recently derived [9,10] for a 
classical Heisenberg spin dimer. The quantum results for arbitrary s are obtained 
using Mathematica® to evaluate the Clebsch-Gordan coefficients. Results have pre- 
viously been obtained for s — \ spin rings of length up to N = 16 by complete 
diagonalization methods [11], and in Ref. [12] some aspects of the high spin limit 
were discussed. 

The present study of the equilibrium autocorrelation function for large values of s 
is timely given the fact that NMR measurements have very recently been performed 
[13] on a dimer molecular magnet composed of Fe 3+ (s = 5/2) ions. Heretofore 
only s = | dimers have been available for NMR studies [14-16]. For comparison 
between theory and experiment it will also be necessary to incorporate molecular 
and single-ion anisotropy terms in the Hamilton operator. This will be the subject 
of a forthcoming article. 



2 The quantum dimer 

The quantum dimer is specified by the Hamilton operator 

H=ih-h = 4v (§ 2 - 1? - li) ; S = li + *2 , (1) 
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where J > describes antiferromagnetic and J < ferromagnetic coupling. 
Throughout this article it is assumed that the spin quantum numbers of both sites 
of the dimer are identical, s 1 = s 2 = s. The eigenstates | S M ) of total spin ~S 2 

S 2 \SM) = h 2 S(S+l)\SM) , S Z \SM) =hM \SM) (2) 

are also eigenstates of the Hamilton operator with eigenvalues E s , which, in the 
absence of a magnetic field, do not depend on the total magnetic quantum number 

M 

H\SM) = - (S(S+1) -2s(s+l)) \SM) = E S \SM). (3) 
2 

Thus the partition function in the canonical ensemble reads 

Z = tr{e -/J ~} = J2( SM \ e_/3 ~ \ SM ) ( 4 ) 

S,M 

= e l3Js(s + l) £ (2iS+1)e -Ms( S+ l) 

and, considering that the Hamilton operator (1) is isotropic, one obtains for the 
unnormalized autocorrelation function 

« 2i(*) • 2i(0) )) = % J2( SM \ s Mt) ■ ^(0) ^ I SM) (5) 

^ S,M 

= 1 eT^~ E s') e-^s \(SM\s lz \S'M')\ 2 . 

Z S,M,S',M> 

The last expressions simplifies when we take into account that only matrix elements 
with M — M' and a difference in total spin not larger than one contribute [17], i.e. 

(SM\ s lz \S'M') = if |S-S'|>1 or M / M' . (6) 

The resulting expression is shown in Eq. (A.l) in the appendix. Note that the an- 
gular frequency spectrum of the spin dimer is given by integer multiples of J/ h, 
which in turn means that the autocorrelation function is periodic in time and that 
the recurrence time r only depends on the coupling J but not on the spin quantum 
number s 

ue^-Sj, S = 0,...,2s T = —j-. (7) 
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Fig. 1. Normalized autocorrelation function Re(C(i)), see Eq. (A.2), for a spin-|-dimer 
for four different temperatures (solid lines). The left panels display our results for the fer- 
romagnetic dimer, the right panels the antiferromagnetic case [20]. The dashed lines show 
the classical result; for details see the next section. 



This is of course true for all Hamilton operators that can be written like the term on 
the r.h.s. of Eq. (1), namely for the spin trimer and the spin tetrahedron. 

Figure 1 shows the autocorrelation function normalized to unity at t — for a 
spin-|-dimer, a system that has been synthesised (Fe dimer) and that is currently 
under investigation [13]. The analytical expression for this autocorrelation func- 
tion is given in Eq. (A.2) in the appendix. One clearly sees that the autocorrelation 
function, which is a superposition of five harmonic oscillations and a constant, is 
dominated at low temperatures by the highest frequency in the ferromagnetic case 
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and by the lowest frequency in the antiferromagnetic case. At higher temperatures 
other frequencies also contribute. One also notices that, independent of tempera- 
ture, the autocorrelation function returns to its initial value after r = 2j ^. 



3 Comparison to the classical dimer 



In order to compare the results of the quantum dimer for different spin quantum 
numbers s with each other and with the classical dimer it is useful to introduce 
normalized spin operators 



In = , = = n=l,2 , (8) 

y^ 2 s(s + 1) 



which depend on s. Note, that 



(9) 



and hence these become commuting operators for s — > oo. Eq. (8) suggests that we 
define a classical Hamilton function H r 



H c = J c e 1 -e 2 , J c = Js(s + l), (10) 

where ei and e 2 are unit vectors (c-numbers). We expect that the thermal proper- 
ties of this classical Heisenberg system will coincide with those of the quantum 
Heisenberg dimer if s 3> 1 except for very low temperatures. This is because the 
spectrum of eigenvalues of e nz is confined within (—1,1) and becomes dense for 
s — > oo and thus coincides with the continuous range of e nz . 

Similarly, if we substitute Eq. (8) in the quantum equations of motion for |i and 
s 2 , we have 



| = —Q |i x e 2 , | = +fi |i x | 2 (11) 
where 



jJs(s + l) J c 

n = —^—t — - = — . . (12) 

h hJs{s + i) 
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This suggests that we prescribe the following equations of motion for the classical 
unit vectors e x and e 2 , 



ei = — Q ei x e 2 , e 2 = +fi ^ x e 2 . (13) 

We emphasize that in these equations f2 is given by Eq. (12). It is expected that 
the autocorrelation function derived using (13) and the canonical ensemble average 
based on H c will coincide in the large s limit with the normalized autocorrelation 
derived from Eq. (A.l). This expectation is in fact confirmed as discussed below. 

Using the fact that the total spin is a constant of motion the classical partition func- 
tion can be derived as [18,19] 



-TjJl dEe ^ E} 

Note that the classical density of states turns out to be a constant in the energy 
interval [— J c , J c ]. This coincides nicely with the quantum density of states which 
can be obtained by counting the discrete eigenvalues per unit energy interval and 
normalizing the density so that its integral gives 1. One can show, starting from 
Eq. (4), that the quantity Z/ (4 s(s + 1)) is in close numerical agreement with Z c 
for temperatures fee T > 0.2 J s(s + 1). This serves to clearly define the classical 
regime for the thermal properties of the dimer. 

For the classical autocorrelation function one finds that [9,10] 



2) j 

sinh(/3 J c ) 



(14) 



C c (t) = - 



1 - coth(/3 J c ) + 



X 



(3J C \ ■ l-exp(-2/3J c ) 

(3J C S 2 ' 



(15) 



J\s S (^-^^j cos (Snt) , 



which can be integrated using error functions of complex arguments, see Ref. [20]. 
In contrast to the quantum autocorrelation function (5), the classical quantity is real. 
The reason is that s^i) • |i(0) is not a hermitian operator for t ^ 0. If one would 

like to construct a hermitian operator, \ • |i(0) + |i(0) • would be 

appropriate. This coincides with the real part of our definition Eqs. (5) and (A.l). It 
is also interesting to note that the imaginary part of ^ gi(t) ■ |i(0) ^ does indeed 
vanish in the high temperature limit. 



In Fig. 2 the dashed curves display the classical autocorrelation function obtained 
from (15) together with the quantum result (solid lines) for three different spin 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

t (2nfi/DI) t (2jrfi/l]l) 

Fig. 2. Normalized autocorrelation function for three different spins at the temperature 
ksT/(J s (s + 1)) = 0.2. The left panels displays the ferromagnetic dimer, the right 
panels the antiferromagnetic one. The solid lines show the quantum result, the dashed lines 
the classical. 



quantum numbers. In order to compare the different correlation functions all spectra 
have been mapped on the same energy interval [— J c , J c ]. Thus the different figures 
show the autocorrelation functions for the same position of the mean excitation 
energy in the spectrum, i.e. the same 



k B T k B T 

"7 — 7 — , 7T = —j— = 0.2 • (16) 

J s (s + 1) J c 

Based on our earlier remark concerning the close numerical agreement of the clas- 
sical and quantum partition functions when k B T > 0.2 J c , we anticipate similar 
agreement for the autocorrelation function in this temperature range. This is con- 
firmed on inspecting the various panels of Fig. 2, which demonstrate nicely that the 
quantum autocorrelation function approaches the classical result with increasing s. 
The most prominent difference between these results is that the classical autocorre- 
lation function does not return to its initial value but approaches a unique non-zero 
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limit, whereas the quantum autocorrelation function is recurrent with a recurrence 
time independent of spin and temperature. This is due to the fact that the classical 
system has a continuous spectrum of excitations in the angular frequency interval 
[0, 2fi] whereas the quantum system possesses a discrete spectrum of excitations 
which are all integer multiples of the lowest one. 
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A The quantum dimer 

Using the matrix properties (6) the unnormalized autocorrelation function (5) can 
be simplified to 



82. 




(A.l) 



2s M=S-1 



+ E E \(sm\s 



S-1M)\ 2 



5=1 M=-S+l 




Taking as an example the case s = 5/2 yields 
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C(t) 



li(*)-li(0) 



(A.2) 



Ii(0)-Ii(0) 

330 + 180 e 5J ? + 84e 9 J/3 + 30e 12 J/3 + 6 e 14 3 p 
+35e J (^ +14/3 ) +35e J ("x +15/3 ) +64e 2J (x +6/3 ) + 64e 2J (-^/3) 

+81e 3j(^ + 3^) +81e 3j(-^ + 4/3) +80e j(Mi +5/3 ) +80e j(-±p +9/ 3) 

+ 55e^ + 55e 5J (-^) 



35 (ll + 9e 5 ^ + 7e 9 J/3 + 5e 12 J/3 + 3e 14 J/3 + e 15 J/3 ) 



Autocorrelation functions for other spin quantum numbers can be evaluated using a 
Mathematica® 3.0 script, that the reader is encouraged to download from our web 
site [20]. 
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